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We study Ivr production in both charged and neutral current neutrino nucleus scattering for 

neutrino energies below 2 GeV. We use a theoretical model for one pion production at the nucleon 

level that we correct for medium effects. The results are incorporated into a cascade program that 

apart from production also includes the pion final state interaction inside the nucleus. Besides, in 

some specific channels coherent n production is also possible and we evaluate its contribution as 

^^, well. Our results for total and differential cross sections are compared with recent data from the 

Cy ' MiniBooNE Collaboration. The model provides an overall acceptable description of data, better for 

NC than for CC channels, although theory is systematically below data. Differential cross sections, 

folded with the full neutrino flux, show that most of the missing pious lie in the forward direction 

^^ ■ and at high energies. 
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I. INTRODUCTION 

A correct understanding of neutrino-nucleus interactions is crucial to minimize systematic uncertainties in neutrino 
oscillation experiments [1| . Most of the new generation of neutrino experiments are exploring neutrino- nuclear scat- 
tering processes at intermediate energies. Recently the MiniBooNE Collaboration has published one pion production 
cross sections on mineral oil by f/i/£u neutrinos with energies below 2 GeV. The data include neutral-current (A^C) 
single 7r° production by iA, and i/^ [2], as well as z^^ induced charged-current (CC) charged pion production Q and 
neutral pion production [j|. These are the first pion production cross sections to be measured since the old bubble 
chamber experiments carried out at Argonne National Laboratory ( ANL) 0, Q and Brookhaven National Laboratory 
(BNL) ,7i]. The latter were measured on deuterium where nuclear effects are small 0, |0j. The main contribution 
to MiniBooNE data comes from ^^C and this poses an extra problem to theoretical calculations as a direct test of 
any fundamental production model is difficult since in-medium modifications of the production mechanisms and final 
state interaction (FSI) effects on the produced pions are important. 

These new data show interesting deviations from the predictions of present theoretical models [lu\ . In Ref. fll| , the 
NuWro Monte Carlo event generator is used to study NC 7r° production data in nuclei. A simple theoretical model 
is used to describe the pion production on the nucleon. No background terms are considered, and the coherent pion 
production is calculated using the Rein-Sehgal model [T^, which is not appropriate at low energies |13l - [l5l |. A fair 
agreement for neutrino induced reactions and a little worse agreement for the antineutrino case are obtained. CC 
single pion production off ^^C for neutrino energies up to 1 GeV has been analyzed in Ref. [16|. The pion production 
theoretical model is here more complete |17| but nuclear effects have been included only in a very simplified manner. 
Only total cross sections were calculated and approximately agree with MiniBooNE data for the 77+ channel and 
deviate more for the tt^ channel. The most comprehensive approach till now, with both a quite complete microscopic 
description of the pion production on the nucleorl]J and the nuclear medium effects, can be found in Ref. }20i] . There, 
the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model is used finding that total cross sections measured by 
MiniBooNE are higher than theoretical ones for neutrino energies above 0.8 ~ 0.9 GeV and obtaining also some 
discrepancies in the differential cross sections. The study of [20] is limited only to the incoherent part of the CC 
induced reaction, and no comparison with the MiniBooNE NC data of Ref. His performed. However, such comparison 
was presented/discussed in NUFACT and NUINT conferences in 2009 [2ll,l2i|. 

In this paper, we address the problem of ;^(P)-induced pion production, both for CC and NC driven processes, 
using a more sophisticated theoretical model, although it is restricted to relatively low pion energies. We start from 
the one-pion production model on nucleons from Refs. [9] and 17] which includes the A resonance mechanisms, but 
also the background terms generated by the leading order chiral Lagrangian. In order to extend the model to higher 
energies above the A resonance region for which it was originally developed, we add a new resonant contribution 
corresponding to the -Di3(1520). According to Ref. [l8|, this resonance, besides the A(1232), is the only one playing a 
significant role for neutrino energy below 2 GeV. Next, we have also incorporated several nuclear medium corrections 
that directly affect the production mechanisms. Apart from Pauli-blocking and Fermi motion we take into account 
the important corrections that stem from A resonance properties modification inside the nuclear medium. Finally, 
the pion FSI is also relevant for the comparison with the experiment, as a number of pions will be absorbed or re- 
scattered, possibly changing their charge, in their way out of the nucleus. For the inclusion of these effects, we shall 
follow Ref. |23| where a simulation code for inclusive pion nucleus reactions was developed. In some specific channels 
coherent tt production in ^^C is also possible and to evaluate its contribution, we use the model of Ref. [1J| but with 
the newer nucleon-to-A form factors extracted in Ref [91], where a simultaneous fit to ANL and BNL data, accounting 
for deuteron effects as well, was carried out. 

The paper in organized as follows: In Sec. |TT] wc introduce our model for tt production induced by viv) at the 
nucleon level, the relevant modifications for in-medium calculations and briefiy describe the implementation of the tt 
FSI. A comparison with MiniBooNE data is shown in Sec. lIIII and the main conclusions of this work are collected in 
Sec. |lVl In Appendix El the details on the D13 contribution to pion production by neutrinos at the nucleon level is 
included. 



^ The model 18] includes the weak excitation of several resonance contributions and their subsequent decay into ttN. The vector couplings 
are taken directly over from the MAID analysis [13 , while the axial couplings are obtained from partial conservation of the axial current 
(PC AC). For the necessary background terms, the vector part is again determined using the MAID analysis as a basis. For the axial part 
(including the vector-axial interference) it was then assumed that it is proportional to the vector part. The proportionality constant is 
adjusted to the old bubble chamber ANL and BNL data, neglecting small deuteron effects. 



II. MODEL FOR ONE PION PRODUCTION INDUCED BY NEUTRINOS 
A. Pion production at the nucleon level 

The starting point is the model of Ref. [l3| for one pion production on the nucleon which is depicted diagrammatically 
in Fig. [TJ It contains the dominant A resonance term and background terms required by chiral symmetry. In total 
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FIG. 1. Model for the W^ N -^ N'-k reaction. We have direct and crossed A(1232)— and nucleon pole terms, contact and 
pion pole contribution, and the pion-in-flight term. We denote these contributions by: AP, CAP, NP, CNP, CT, PP and 
PF, respectively. 

we have direct and crossed A(1232)— (first row) and nucleon— (second row) pole terms, contact and pion pole 
contribution (third row) and finally the pion-in-flight term. The background terms are the leading contributions 
of a SU(2) nonlinear a model. Those were supplemented with well known form factors in a way that respected 
both conservation of the vector current (CVC) and PCAC hypotheses. Their contribution is sizeable even at the 
A(1232)resonance peak and it turns out to be dominant near pion threshold. 

For the nucleon to A weak transition matrix element the form factor parametrization of Refs. [2J, [2^ was 
takero in [17]. A total of four vector and four axial form factors are needed. The three vector form factors 
C^(g^), CY {q"^) , C^ {q'^) were determined from photo and electroproduction data w;hereas C^(g^) = from CVC. The 
vector form factors from Ref. [2a] were used. For the axial part, Adler's model [2,7] in which C^{q^) = 0, Cf{q^) = 
— C^(q^)/4, and PCAC, that requires C^iq"^) = —C^{q^)M'^/{m'^ — q^) with M, m^r the nucleon and pion masses, 
were used. Thus, there is only one axial form factor to be determined, namely the dominant C^{q^). For the latter, 



the parametrization of Ref. [28| was assumed in [17[ and the unknown parameters C^ (0) and M^ia were fitted to the 



fiux averaged, W < 1.4 GeV, i^^p — > fi^pir^ q^ differential cross section measured at ANL [6f. Here, W stands for the 
final pion- nucleon invariant mass. 

While the results for different total and differential cross sections were in good agreement with ANL data [a, |g] , the 
total cross sections were smaller than the experimental data measured at BNL [3]. BNL and ANL data seemed to be 
incompatible. It was pointed out in Ref. [29] that this problem might originate from two factors: First, both ANL and 
BNL data were measured on deuterium. Deuteron structure effects in the j/^d — > fi^A'^'^n reaction were estimated 
in Ref ^ to produce a reduction in the cross section from 5 — 10% and most analysis, including Ref. [fj], neglected 
that effect. Second, both experiments suffered from neutrino flux uncertainties that in Ref. 29] were estimated to 
be 20% for ANL and 10% for BNL. Following the work of Ref. [3, a combined fit of C^(q^) to both ANL and 
BNL pn'^ data, including in both cases full deuteron effects and the flux normalization uncertainties, was carried out 
in [g]. In this latter fit, a simpler pure dipole parameterization C^{q^) — C^(0)/(1 — q^ /M\^'^ was used obtaining 



Note that the Cg sign is quoted incorrectly in Ref. [24| (see comment in Ref. |24l . |25| ) . 
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FIG. 2. Comparison of the theoretical model Q] results (solid line) to ANL 6j and BNL [7[ experimental data. Theoretical 
68% confidence level bands are also shown. Data include a systematic error (20% for ANL and 10% for BNL data) that has 
been added in quadratures to the statistical published errors. The theoretical results and ANL data include a. W < 1.4 GeV 
cut in the final nN invariant mass. 
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FIG. 3. I?i3(1520) contributions to W+N -> N'n. We have direct (DP) and crossed (CDP) Di3(1520)-pole contributions. 

C^(0) = 1.00 ± 0.11 and Maa = 0.93 ± 0.07 GeV. The results from that fit, compared to ANL and BNL data, are 
shown in Fig. [2] The above mentioned parameterization for C^{q^) is the one we are using in the present calculation. 

For the present calculation, and in order to better compare with MiniBooNE data, we need to extend the model 
up to 2 GeV neutrino energies, where higher mass resonances could play a role. According to Ref. [l^, and apart 
from the A, only the I?i3(1520) resonance gives a significant contribution in that region. We shall include in our 
model the two new contributions depicted in Fig. |31 All the details of the calculation of these diagrams are given in 
Appendix]^ We would like to remark that as the D13 has isospin 1/2, it does not contribute in the p7r+ channel and 
thus, it does not affect the previous fit. 

The differential neutrino-nucleon cross section with respect to the pion energy Et, and the angle between the pion 
momentum and the neutrino beam 0^ is given for a CC process by |17| 
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(1) 
Gp is the Fermi decay constant, k, k' and fc^r are the neutrino, final charged lepton and final pion four-momenta respec- 
tively. Besides, q — k — k' is the four-momentum transferred, and En = \J M"^ -I- p^ , E^^i — d A'P -\- (pV + q — k^^Y 
are the initial and final nucleon energies. The lepton and hadronic tensors are given by 

•'-/^o-(K, k ) ^ "'/j"'(T I ^a-l^fj, k ■ K g^^ + I'^fj.aa/S'^ "^ (2j 

W^'^{pN,q,k^) =^ {N'7r\ji!^40)\N) (iV'^b-+(0)|iV)* (3) 



with eoi23 = +1 and the metric g^i, = (+, — , — , — ), and where in the hadronic tensor a sum over final spins and an 
average over initial ones is done. The Jcc+i^) current contains all contributions in Fig.[Tl see Refs. [y,[i3] for details, 
plus the new D13 contributions of Fig. |3] that are given in App endix El A similar expression is obtained for an NC 
process in which we have a neutrino in the final state (see [17|). 



B. In- medium production 

For incoherent production on a nucleus we have to sum the contribution to the cross section of all nucleons in the 
nucleufl Assuming the nucleus can be described by its density profile p{r) — Pp{r) + Pn(r), and using the local 
density approximation, the initial differential cross section at the nucleus level for a pion production channel N'lr^ 
prior to any pion FSI, is then 



N—n.p 



fd'^ E 2 / ^ OiE^ir) - E.) OiE. + ,° - i.. - Ef (.)) '^^f ^ ^f "^ , (4) 



where Ep{r) — \ M"^ + {kp{r)Y ^ being kp{r) = (37r^p7v(^))^^'^ the local Fermi momentum for nucleons of type N . 
To compare with experiment, we have to convolute the above expression with the neutrino flux $(|A:|) 

-4-- = /d|fc>(ifci)4. /.^.. 5:2 Z;^ .(<(,) ^ 

acosUTT cIEt^ J J ^f-^ J (27r)'^ acosu-^dh-n 
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(5) 
From there, we obtain 



d\k\^i^r^drdco^B^dE^ Jl^^^J^^^f a'; V a' y f K )) dcosO^dE^ 



(6) 



Apart from modifications discussed in what follows, the above differential cross section is used in our simulation code 
to generate, in a given point inside the nucleus and by neutrinos of a given energy, pions with a certain charge, energy 
and momentum direction. 

Defining P = q — k^ (the four momentum transferred to the nucleus) and writing d^pN = d cos i}n dipN \pN\ENdEN, 
where the angles are referred to a system in which the Z axis is along P, we can integrate in the i^at variable using 
the energy delta function present in '^^'^^^ e dE ^ "^^^ final result is 



d\k\ATTr^ dr d cos 



'^^ = $(1^1) [ dn'dE'ik'il V ,f/, i^-L g«(r)-g 



)e{-p')0{p'>)c,,,{k,k') 



2n i-E'^ir) 

d^N / dENW^'^PN,q,K) 



COS l^iv— COS 1^9^ 



(7) 



where 






To speed up the computational time, we approximate the last two integrals in Eq.([7l) by 
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(9) 



where pN is evaluated at the value Ej^ = {Ep (r) +£)/2, (middle of the integration interval), with the corresponding 



cos-!^^ (that deduced from Eq. ^ using E^ = Em and \pn\ = J Ef^ — M"^ ), and ^at is set to zero. Similar 
approximations were done, and shown to be sufficiently accurate, in Refs. |3CH33l| to study total inclusive and total 
inclusive pion production in photon and electron nuclear reactions. In the study in Ref. [3J] of total inclusive neutrino 
induced cross section this kind of simplification was also used. We have checked that the approximation of Eq. (|9| 



For coherent production one should sum amplitudes [13 



induces uncertainties at most of 5%, independently of ^jy. Other choices to fix (^at produce small variations of the 
order of 1-2%. With this approximation, we find 

da 



d\k\ATTr^drd cos 9,rdE^ J ^ n~^,^^^^ \P \\>^\ 

xCf,„{k,k')W"'{pN,q,kA. (10) 



Eq.lfTU]) includes explicitly Fermi motion of the initial nucleon and Pauli blocking of the final nucleon but there 
are other important in-medium corrections, that we discuss in the following, that have to be included. The above 
expression is equivalent to Eqs. (2) and (25) of Ref. [3J|, where among others, the IplhlTr excitations contribution 
to the total inclusive neutrino-nucleus cross sections was evaluated. 

C. In medium corrections to pion production and final state interaction 

Given the dominant role played by the AP contribution and since A properties are strongly modified in the nuclear 
medium |13ll32l. [35l - l4]| a more proper treatment of the A contribution is needed. Here, we follow Ref. [S^l and modify 
the A propagator in the AP term as 

1 1 1 



pI-MI + iMATA Vs + Ma y^ - Ma + i(rP^"'V2 - LtiSa) ' 

with s = p\, pPauii ^jig fi-gg /^ width corrected by Pauli blocking of the final nucleon, for which we take the expression 
in Eq.(15) of Ref. [42|, and IuiSa the imaginary part of the A self-energy in the medium. For the mass we shall keep 
its free value. While there are some corrections to the mass coming both from the real part of the self-energy and RPA 
sums, together they induce changes smaller than the precision in the present experiments and the uncertainties due 
to our limited knowledge of the nucleon to A transition form factor C^{q^), see discussion in Sec. II. E of Ref. [34] . 
The evaluation of Ea is done in Ref. [36| where the imaginary part is parametrized as 

-I,nE..C«(^)\c..(iL)\a..(£' 

with po = 0.17fm~'^. The CQ,a, Ca2,P and Cas,^ coefficients can be found in Eq.(4.5) and Table 2 of Ref. (SQ- 
They are parametrized as a function of the kinetic energy of a pion that would excite a A of the corresponding 
invariant mass and are valid in the range 85MeV < T^^ < 315 MeV. Below 85MeV the contributions from Cq and 
Ca^ are rather small and we take them from Ref. '42], where the model was extended to low energies. The term with 
Ca2 shows a very mild energy dependence and we still use the parameterization from Ref. [36i] even at low energies. 
For T^ above 315 MeV, we have kept these self-energy terms constant and equal to their values at the bound. The 
uncertainties in these pieces are not very relevant there because the A — >■ Nt: decay becomes very large and dominant. 
The terms in Ca2 and Ca^ are related to the two-body absorption WNN -^ NN and three-body absorption 
WNNN — > NNN channels respectively. On the other hand the Cq term gives rise to a new WN — ^ Ntt contribution 
inside the nuclear medium and thus it has to be taken into account beyond its role in modifying the A propagator. 
This new contribution has to be added incoherently and we implement it in a approximate way by taking as amplitude 
square for this process the amplitude square of the AP contribution multiplied by 

CQJP/Por ,^^^ 

pFroc/2 • ^ '' 

Our final model for production is then given in Eg. ljlOp with the modifications in the hadronic tensor just mentioned. 

Once the pion is produced inside the nucleus, it starts propagating and it suffers interactions with the medium. To 
evaluate them we follow Ref. [23| , where a computer simulation code was developed to describe inclusive pion nucleus 
reactions (quasielastic, single charge exchange, double charge exchange and absorption). We take into account P- and 
S'-wave pion absorption, and P-wave quasielastic scattering on a single nucleon. The P- wave interaction is mediated 
by the A resonance excitation where the different contributions to the imaginary part of its self-energy give rise to 
pion two- and three-nucleon absorption and quasielastic processes. After a quasielastic interaction the pions change 
direction and may change charge. The intrinsic probabilities for each of the above mentioned reactions were evaluated 
microscopically as a function of the density and we use the local density approximation to evaluate them in finite 



nuclei. In between collisions the pions are treated as classical particles and in the present calculation we shall assume 
they propagate in straight lines. All details of the simulation can be found in Ref. ^2^]. We should remark that this 
approach to tt FSI has been extensively and successfully used in the study of many processes such as hypernuclear 
decays 1431. tt absorption 44 [, photon [31, 45] and electron fSS"] induced nuclear processes (like tt production) or muon 



capture |46[. Furthermore, it has been extensively used in their analysis by the successive Kamiokande collaborations. 
See, e.g. Refs. ^^ 



D. Coherent production 

In some of the channels there is a contribution from coherent tt production. We evaluate it with the model of 
Ref. jl^, using the C^{q'^) form factor obtained in Ref [9|. Cross sections for the T2K and MiniBooNE fluxes 
obtained with this new form factor were given in [4l| . The model for coherent pion production is based on the 
microscopic model for pion production off the nucleon of Ref. [T^j that as already mentioned, besides the dominant 
A pole contribution, takes into account the effect of background terms required by chiral symmetry. Our coherent 
production model does not incorporate the -D13 contribution that we expect to produce only a small correction. 

The main nuclear effects, namely, medium corrections on the A propagator and the final pion distortion, are 
included. As found in similar calculations [33, llo, ImI , the modification of the A self-energy inside the nuclear medium 
strongly reduces the cross section, while the final pion distortion mainly shifts the peak position to lower pion energies. 

We should stress that the model of Ref. [ij] is more reliable than the Rein-Sehgal approach [l^] for the energies 
of interest in this work. In particular, it greatly improves on the description of angular distributions of the outgoing 
pion with respect to the direction of the incoming neutrino [13l |l5| . 



III. RESULTS FOR PION PRODUCTION IN NEUTRINO NUCLEUS SCATTERING AND 

COMPARISON WITH MINIBOONE DATA 

In this section we compare our predictions with data recently obtained by the MiniBooNE Collaboration for i/^ 
induced CC [Sl, l^] and v^i/v^^ induced NC Q] pion production in mineral oil {CH2). 



A. CC production 

We start by showing results for the total unfolded cross sections for charged current one-pion production. In Fig. 21 
we compare our results with the data by the MiniBooNE Collaboration for a final 7r+. We take into account the 
contribution on ^^C and that on the two hydrogens. There is also a small coherent contribution on ^^C that we have 
evaluated as described above. Our total result agrees well with data at low neutrino energies, but is below data for 
neutrino energies above 0.9 GeV. 

The D13 contribution is only noticeable above E^ — 1.2 GeV. At around Ei, — 2 GeV, it makes some 8% of the 
total contribution. It will play a minor role for observables that are convoluted over the MiniBooNE neutrino flux 
as the latter peaks at around 0.6 GeV. On the other hand, in the whole energy range shown, the Cq term (Eq. pTjl ) 
produces changes always smaller than 10%. In Fig.[5lwe show the effect of varying C^(0) within the uncertainties in 
its determination, using BNL and ANL data, in Ref. ^. We find effects at the 10% level. The highest value seems to 
be favored by the MiniBooNE data. We should remind here that if fiux uncertainties were ignored, BNL and ANL 
data would not be compatible [29]. Taken at face value, the figure suggests that MiniBooNE data would sup por t 
BNL results. This would also imply larger values of C^(0) closer to the PC AC prediction, as suggested in Ref. f20| . 
However, one should be cautious given the number of uncertainties in both experiment and theoretical models. 

Similar results, consistently below data, are obtained for the case of CC tt" production, see Fig. |6l As in the 
previous case, the role of the Di^ resonance or the Cq corrections are small. 

In Fig. [71 we compare the differential -^jjr- cross section for CC l7r+ production by f^, calculated with the MiniBooNE 
flux from Ref. [3|. In the left panel, we show the different contributions to the full model, including the coherent 
part. The model predicts less high energy pions than the experiment for T^r above 0.15 GeV. The combined effect of 
quasielastic scattering and pion absorption through A excitation depletes the 0.15 ^ 0.4 GeV region. As a result the 
strength moves down to 0.8 GeV. The shape of the calculated cross section is dominated by the ^'^C contribution and 
both the coherent part and the hydrogen contribution peak at higher energies but, they are too small to compensate 
for the missing high energy pions. On the right panel we also show the effect of not including the Cq correction of 



u 

N 

B 

u 

'o 

D 



2H 
Coh 

'^C+2H+Coh 

'^C (no C„)+2H+Coh 



V N ^ U. 71 N' 



C (no N*)+2H+Coh 
MiniBooNE 



140 
120 
100 

80 

60 

40 

20 


U.3 J 

E^ (GeV) 

FIG. 4. Ivr"'" total production cross section for Ufj, CC interaction in mineral oil. Short-dashed line: ^^C contribution. Dotted 
line: H2 contribution. Double-dashed dotted line: Coherent contribution (see main text). Solid line: Total contribution. 
Long-dashed line: Same as solid line but without the Cq contribution of Eq. (|lip . Dashed-dotted line: Same as solid line but 
without the contribution from the D13 resonance contribution. Data taken from Ref. [3|. 
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FIG. 5. Itt"*" total production cross section for u^ CC interaction in mineral oil. Solid line: Our full model with C5 (0) = 1. 
Dashed line: Full model with C^(0) = 1.1. Dashed-dotted line: Full model with C^(0) = 0.9. 



Eq. (fTTj) . We can see, that it little affects the cross section and/or the shape of the pion kinetic energy differential 
distribution. 

From the above discussion we might expect a better agreement with data when FSl effects are neglected. Indeed, 
as can be seen on the right panel of Fig. [3 the experimental shape could be reproduced by artificially removing the 
FSl of the pion. Of course, this has little sense as FSl is really there. Removing FSl completely would lead to total 
cross section values that are too high at low neutrino energies, where the theoretical model is more reliable. FSl 
effects could be reduced by considering the so called "formation zone" [ll[, that among other effects includes the 
propagation of the A before decaying into a ttN pair. Of course, the "formation zone" could be adjusted to reproduce 
data. These kind of modifications of the FSl could be difficult to justify, they might be in confiict with much other 
phenomenology and may somehow serve to hide our ignorance on the relevant dynamics. Even when they could help 
reproducing some observable, if we lack a correct understanding of the physical mechanisms responsible for them, 
they might lead to wrong predictions for other observables sensible to other kinematics, dynamical mechanisms or 
nuclear corrections. 
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FIG. 6. Itt" total production cross section for i^^j CC interaction in mineral oil. In this case there is no contribution from the 
iif2 and there is no coherent contribution. Captions as in Fig. [l] Data from Ref. [J]. 
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dotted line in the left panel) contributions. Data from Ref. [3]. Note that the coherent (left panel) and D\z resonance (right 
panel) contributions are very small. The curves obtained when these effects are neglected can hardly be distinguished from the 
red solid line that stands for the total contribution. 



In FiglHl we show the flux-folded differential 4^ and 



cross sections for CC Itt" production by Vn that we 



"COS Wtt ^ J ii- 

compare to data. For that we use the neutrino flux reported in Ref. [J], that extends from 2 GeV down to 0.5 GeV 
neutrino energies. In this case, there is neither contribution from the hydrogen nuclei nor from coherent production. 
Once more, the FSI effects are clearly visible in both distributions and their artiflcial exclusion leads to a better 
description of the high momentum tail of the j^L distribution. Because of the FSI some pious are absorbed, but other 
ones are scattered and loose to nucleons part of their energy. FSI is essential to fill the low momentum part of the 
distribution. From the angular distribution, we see that the pion production off the nucleon model of Ref. [13| leads 
to a forward peaked cross section. FSI, through quasielastic collisions tends to soften the curve and leads to a better 
description of the backward scattering. From our point of view, these figures might indicate that some mechanism 
for pion production, that provides forward high energy pious could be missing in our theoretical scheme. 
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FIG. 



Flux-folded differential -j^ (left panel) and 



oil. Captions as in Fig. [71 Data from Ref. [J]. 



(right panel) cross section for CC Itt" production by u^ in mineral 



B. NC production 

In these channels, we compare our model results to data from Ref. [2|- In each case we use the different neu- 
trino/antineutrino fluxes reported by the MiniBooNE collaboration. In Fig. ^ we show the j^ differential cross 

section and the different contributions coming from ^^C, H2 and the coherent production on ^^C. Both for neutrino 
and antineutrino reactions, we see that the model agrees better with data than in the CC case. We still obtain cross 
sections below data in the 0.25 ~ 0.5 GeV/c momentum region. We also see that the role of the D13 is negligible 
while the effects of the Cq term are small. 
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FIG. 9. Flux-folded differential -j^ cross section per nucleon for NC Itt*^ production by ly^ (left panel) and i/^ (right panel) 

in mineral oil. Short-dashed line: C contribution. Dotted line: H2 contribution. Double-dashed dotted line: Coherent 
contribution. Solid: Full model result. Long-dashed line: ^^C contribution without the Cq term (Eq. HI}). Dashed-dotted 
line:^^C contribution without the D13 term. Data from Ref. [2|. 



In Fig. [To] we show now the flux- folded differential 



dcr 



cross section. The full model agrees better with data in 



the antineutrino case where our results are within error bars except in the very forward direction. The role and size 
of the coherent piece is crucial for the agreement. 

There is some deficit for cosOt^ > for the reaction with neutrinos but the agreement, as it was the case for the 
■j^ differential cross section, is better than in the corresponding CC reaction. In fact, a minimal enhancement of the 
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panel) in mineral oil. Captions as in Fig. |9l Data from Ref. 0]. 



cross section per nucleon for NC Ivr" production by u^ (left panel) and i^fj, (right 



coherent process and/or the hydrogen contribution, that could be obtained by a larger value of C^ (0) also suggested 
by the CC results, could lead to a better agreement. 



IV. CONCLUSIONS 

We have extended the model for pion production by neutrinos on nucleons from Refs. |9l. llTJ] by the inclusion of a 
new term related to the production and decay of the L)i3(1520) resonance. The resulting model is expected to give a 
fair reproduction of experimental data on nucleons for neutrino energies up to 2 GeV, always with the uncertainties 
associated with the poor knowledge of some relevant form factors. 

Including nuclear medium corrections that affect both the elementary tt production mechanisms and the pion FSI, 
we have studied pion production by neutrinos in mineral oil in order to compare with recent experimental total 
and differential cross section measurements. The model provides an overall acceptable description of data, better 
for NC than for CC channels. In the CC channels, the predicted total cross sections are below data for neutrino 
energies above 0.8 ~ 0.9 GeV. This result is in agreement with other theoretical calculations [lg,[20]. Differential cross 
sections, folded with the full neutrino flux, show that most of the missing pions lie on the forward direction and at 
high energies. This might suggest the need of further production mechanisms. NC channels show a better agreement, 
although theory is also below data. We find that the role of the coherent tt production is essential to properly describe 
the angular distributions in the NC channels. We also find that flux unfolded results seem to support a large value 
of C^(0), closer to the PCAC prediction than the value obtained from previous combined fits to the ANL and BNL 
bubble chamber data. Actually, as already discussed in Ref. ^J, it seems that MiniBooNE data would be better 
described using a C^(0) value fitted to BNL data alone. 
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Appendix A: D13 resonance contribution to neutrino-nucleon scattering 



In this appendix we give full details of the -D13P and CDi^P hadronic matrix elements of the Jcc+{0) weak current, 
depicted in Fig. |31 Those we shall call for short j^^_^_\dp and j^c+Icdp- 
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The matrix element for the direct (DP) contribution is given by 



Ji^c+lop^iC god ^ cos 



k2 



'pJj-MJj + iMdTd 



(Al) 



pn^ 



with Mo = 1520 MeV the mass of the I?i3 resonance. For go we take go — 20 GeV~^ which resuhs from a fit of the 
Di3 — > Ntt decay width. For the latter we take 61% of 115 MeV. The width has two main contributions, the Ntt and 
the Att channels. In the propagator we shall use 



i£)-l_D +'^Dy 



where for F^'^ we take 



with s = p^ and Ip^rl — 2 r '™" ' being A(a, b, c) ^ a? + b"^ -\- c^ — 2ab — 2ac — 26c. 

For F^'^ we assume an 5— wave decay and take 

\p:\ 



r^"" = 0.39 X 115 MeV ,^_ 61(^5 - M - m^) 



\P^' 



with IpJI = 
Besides, 






and \p;°- 



2Mn 



PUpd) = -i^D + Mo) 



ga/3 - 37.7/3 - 3 ^^2 



' 3 M. 



D 



and 



T^^{p,q) 



^ {g^'i ~ q'r) + % (/^<z • PD - q'p'o) + % (g'^q ■ p - q'p') + cW 



% h'H - q'r) + % (g'^q ■ PD - q^'o) + C^'' + S^'?'' 



M2 



M2 



The axial form factors are taken from Ref. [^ 

-2.1 



1 



O3 -O4 -U, O5 - (1 _ ^2/^2 )2 1_ ^2/(3^2)^ 



Al„2\ 



, Ci{q')^Ct{q') 



M2 



9 9 ' 



75, PD=P + q- 



Ma = 1 GeV, 



M7 -^ y / VJ^^A 

while for the vector ones we fitted the form factor results in Ref. [54| to get 

-2.98 ~v 4.21/1?^ ~y -3.13/Dv 



a 



^ [l-g2/(1.4M2)]2' ^4 1 _„2/^^7A-f2^' C's 



l-g2/(3.7Af2)^ 



— r^ = 

2//n /ioA/r2\' ^6 



1 - g2/(o.42M2 ) ' 



with Mv = 0.84 GeV and Dy = (1 - q^/M^f. 
For the crossed CDP contribution we have 



3rr+\CDP 



V 3 pj, - M|) + iMoTo 



t'\rf^0 iJ 



u{p')T'i^^ip',~q)P/^^{poh,uip) 



with 



and 



PD =p - q , C 



CDP 



>^ /^' 



Orr/^M t^i 



r%'{p',-q) = i'n'ip'^-<iWi 



1 p-K^ 

m:^ 



lt^,o 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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From isospin symmetry we further have [ij 



{p7r-\C40)\p)^{n7r+\j^Jn), 



("Tr iJc 



i) = {pTT+\j^Ap). 



For NC neutrino and antineutrino induced reactions the contribution from the D13 has isovector plus isoscalar 
parts. The isovector (IV) part is related to the CC processes and is given by (see discussion on Sec. Ill of Ref. [17]) 



(P»°l«c,Iv(»)IP) = 



'OS Hr^ '^ 



\/2 COS 6c 



l-2sm^ew) 



{p7r+\V^,A0)\p) + {n7r+\V^:A0)\n) 



(p^+|A^^,(0)b) + (n^+|A^^+(0)|n)]} 



(^'^+b;':c.iv(0)b) = -:3^{ (1 -2sm'e^)[{pn+\V^^,mp) {nn+lV^^^, 



cos 6c 



{pn+\A'^^^,{0)\p) - {nn+\A^:^AO)\n)]} 



{P^-\jLwmn)=-{nn+\j^^^^,^mp). 

where we have written j^ + (0) — V^^ (0) — A^ + (0) being F''+ (0) and A'^, + (0) respectively the vector and axial part 
of the current. 

The isoscalar current is given in terms of the isoscalar part of the electromagnetic current as 

J::c.is(O) = -4sin2 0H^s,^^jg(O), 

and one has (see discussion on Sec. Ill of Ref. |17| ) 

{nn+\s':^jsmp) ^ (P7r-|4„^js(0)|n) = ^/2 (p7r"|4„_j3(0)|p) = -^2 (n^°|sL_js(0)|n). 

The direct contribution to ^[('T-Tr^ls^jjj jg(0)|ri) — (ptt^IsJ^j^^ jg(0)|p)] is given by 



2L\""' l''om,ISV";i'V VP" I'^cmJS 

1 k. 



-igo- 



y^pl^All+iMDTD 



u{p'hbPafi{PD)Tji ^^{p,q)u{p), pD=p + q, 



and 



with 



D 



a 



VJS 



a 



y,is 



a 



l^.IS 



M 



b'H - i^r) + ^ (g'^i ■ PD - q^p'n) + ^ (g'^i ■ p - q'pn + cV'g^' 



^v.is -1-21 ^ 

3 [1-92/(1-4M2)]2' ^^ 



IVP 



^yjs 0M5/D^ 



AP 



1 - gV(3.7M^) 



, C; 



y,is 



0.395/1?^ 



1 - gV(0.42M^) 



2\ ' ^6 



cP' = 0, 



with their values at g^ = extracted from Ref. [52]. 

The crossed contribution to ^[{nn'^\s'^^^g{0)\n) — (p7r°|s(^jj^ jg(0)|p)] is given by 



k^ 



^go 



y/Spl-Ml+iMoTo 



u{p')r'/^^^ {p\ -q) P;i>JpDh,u{p), PD=P'- q, 



with 



-Ay.ISufl/ / \ Oi-nV.lS 13 u/ I Mt 
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